Effects of impurities and vortices on the low-energy spin excitations in high-T c materials 
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Abstract 

We review a theoretical scenario for the origin of the spin-glass phase of underdoped cuprate materials. In particular it is shown 
how disorder in a correlated of-wave superconductor generates a magnetic phase by inducing local droplets of antiferromagnetic 
order which eventually merge and form a quasi-long range ordered state. When correlations are sufficiently strong, disorder is 
unimportant for the generation of static magnetism but plays an additional role of pinning disordered stripe configurations. We 
calculate the spin excitations in a disordered spin-density wave phase, and show how disorder and/or applied magnetic fields lead to 
a slowing down of the dynamical spin fluctuations in agreement with neutron scattering and muon spin rotation (/iSR) experiments. 
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1. Introduction 
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From the phase diagram of the cuprate superconductors it is 
evident that as a function of doping, antiferromagnetic (AF) or- 
der gives way to li-wave superconductivity (dSC). This has led 
to the natural view that these two ordered phases generally com- 
pete in these materials. At an overall global scale this certainly 
seems to be the case. However, when including disorder the op- 
posite can also occur: dSC can function as a catalyst for induced 
spin-density wave (SDW) order| 1]. This takes place because of 
the generation of low-energy impurity resonance states which 
are spin split by local magnetic ordering ill. A qualita- 
tively similar effect happens in the vortex phase generated by 
an external applied magnetic fieldjl]. The of-wave symmetry of 
the superconducting state is crucial for this cooperative effect 
to occur between SDW and dSC phases. 

Experimentally, it is well known that disorder may cause 
slowing down and eventual freezing of spin fluctuations @, H 
0,0]. For example, substitution of nonmagnetic Zn ions for 
Cu in near-optimally doped La2- A Sr A Cu04 (LSCO) has been 
shown to shift spectral weight into the spin gap, and eventually, 
for enough Zn (~ 2%), generate elastic magnetic peaks in the 
neutron response |7|. In YBa2(Cui_ v Zn y )306.97 with y = 2% a 
similar in-gap Zn-induced spin mode was observed||9|. Upon 
increased temperature T the elastic signal decreases and even- 
tually vanishes near T c which is similar to an equivalent disor- 
der signal in Zn-free LSC0[6, 10], These results of disorder- 
induced freezing of spin fluctuations are consistent with /jSR 
data on underdoped cupi-atesE! E2 H H H • 

Recent neu- 
tron scattering off detwinned YBa 2 Cu 3 6 .6 (YBCO) with 2% 



Zn found induced short-range magnetic order and a redistri- 
bution of spectral weight from the resonance peak to uniaxial 
incommensurate (IC) spin fluctuations at lower energies lflal . 

Applied magnetic fields introduce vortices into the system 
and cause much of the same phenomenology as described above 
for the disorder. Theoretically this can be ascribed to the fact 
that vortices in dSC also generate low-energy resonance states 
which may favorably split due to electronic correlations and 
generate local SDW orderfllHOdSIIIllzl. For sam- 
ples without static SDW order, experiments show that at low 
magnetic fields the vortices slow down the spin fluctuations 
and generate an in-gap modelSlliQ. For larger fields the 
ground state exhibits static local SDW order which disappears 
above T c . Enhancement of static order was reported first in neu- 
tron diffraction experiments[6, 26] on LSCO, with a correlation 
length of several hundred A, and has been confirmed in other 
underdoped cuprates 11271 1281 129 . 30]. 

The freezing of spin fluctuations is also relevant even in the 
absence of substitutional disorder or external applied magnetic 
fields, as evidenced by the presence of a spin-glass phase in 
the underdoped regime. This seems to apply to both "clean' 
cuprates like YBCO where quasi-static SDW order is found| 3 1 
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33J,[34J], and to intrinsically disordered materials like LSCO 
where the static spin correlations are long-range, and persist 
for a large doping range well into the dSC dome ll3~5[l36ll . The 
size of the spin-glass phase in temperature and doping is clearly 
enhanced by disorder. 

Here, we review a picture that has emerged for the de- 
scription of the induced magnetism in underdoped cuprate 
materialstlll 14. l37l 138. 39, 40]. We focus on the disorder- and 
field-induced SDW phase coexisting with dSC, and calculate 
the associated dynamical spin susceptibility. 
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2. Model 



3. Results and Discussion 



The basis for our model analysis is the BCS Hamiltonian for 
a dSC with orbital coupling to an applied magnetic field B, to 
which we add site-centered disorder and a local Hubbard repul- 
sion treated in an unrestricted Hartree-Fock approximation: 
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Here, c. creates an electron on site i with spin a -\, | on a 
N x N square lattice. The hopping matrix elements between 
nearest and next-nearest neighbor sites are denoted by fy = t 
and tjj = t' , respectively. The magnetic field is included through 
standard Peierls phase factors (pn = (tt/Oo) C' A(r) ■ dr, where 
<E>() = hc/(2e) and A(r) = B(0, x) is the vector potential in Lan- 
dau gauge. The chemical potential fi is adjusted to fix the elec- 
tron density n — YiA n i) = 1 — JC, where x is the hole concen- 
tration. The <f-wave pairing amplitude Ay is determined by the 
strength of an attractive nearest-neighbor interaction Vd- The 
impurity potential V"" p arises from a concentration n„, v of non- 
magnetic pointlike scatterers at random positions. All fields, 
Ay, the local charge density (n,), and the local magnetization 
(cr|) are calculated self-consistently from the Bogoliubov-de 
Gennes (BdG) equations!! . 21 ] . In the following we fix pa- 
rameters to ?' = -0.4?, V d = 1.34?. 

The model ([TJ has been used extensively to study the com- 
petition between dSC and SD W p hases[41, the origin of 
field-induced magnetizationll ll 1 1 7l 1 1 9l \2(X |43 l| . and moment for- 
mation around nonmagnetic impuritie s ikl 1 3 7l 1 3 9ll . In the case of 
many impurities, Eq.(Q]) was used to explain the origin of static 
disorder-induced magnetic droplet phases|37, 39, 44, 45], and 
to study how these may form a quasi-long range ordered SDW 
phase. The underlying idea is that when magnetic droplets 
begin to overlap, either when the magnetic correlation length 
grows at low doping, or when the concentration of scatterers 
is increased, the effective interaction between droplets allows 
them to align their staggered patterns coherently 1371 15311 . More 
recently Eq.([TJ was used to obtain semi-quantitative description 
of the ^-dependence of the elastic neutron response in LSCO 
in an applied magnetic fieldQ],!^]. Finally, Eq.([T]i has been used 
to explain various transport measurements ll39l 14611 such as, for 
example, the experimental observation of a non-universal low- 
T limit of the thermal conductivity k{T) 11471 14811 in terms of 
a disorder-induced SDW phasel39J]. The SDW droplet phase 
has an enhanced scattering rate while maintaining the same 
low energy density of states (DOS) as the homogeneous case, 
thus breaking the cancellation between the residual quasipar- 
ticle DOS and relaxation rate which gave rise to the universal 
k(T) in the first place 149115011 . Thus, for underdoped cuprates it 
is not k(T) but rather the (spatially averaged) low-energy DOS 
which is universal, in agreement with STM experiments! 51]. 
Theoretically a universal low-energy DOS can be traced to sup- 
pressed c harge m odulations caused by the electronic Coulomb 
repulsioniJIMIil. 



The Hamiltonian (Q]) supports both a correlation- and 
disorder-induced SDW phase. Specifically, in the clean case 
above a critical repulsion U c i a global stripe phase is the fa- 
vorable state, and disorder acts mainly to scramble the stripes. 
Below U C 2 the ground state is a homogeneous dSC but non- 
magnetic disorder or vortices may locally induce SDW order if 
U > U c i where U c \ is another critical interaction streng thflli. 
In this low-f/ regime, both disorder and the of- wave symmetry 
of the pairing are crucial for generating static magnetism. 

In the absence of a magnetic field, Fig. Q~f a - C ) show the mag- 
netization in real-space arising from a weak disorder potential 
from the dopant ions l37ll . Figure [TJa) shows that not all impu- 
rities in the correlated system need to "magnetize" for a given 
U; in the disordered system, the effective criterion to drive the 
impurity through the local magnetic phase transition is different 
for each impurity. Increasing the repulsion U then increases the 
concentration of impurities which induce a local magnetization 
droplet, as shown in Fig. |TJb) andQJc). With further increase of 
U, the system evolves from a state with dilute non-overlapping 
AF droplets to connected spin textures. A similar process takes 
place in ID where Shender and Kivelsonj53|] pointed out that 
the interactions between impurities in a quantum spin chain are 
non-frustrating: if an impurity creates a local AF droplet, a sec- 
ond one can always orient itself to avoid losing exchange en- 
ergy. In 2D this continues to apply for spin models with nearest 
neighbor exchange, but may break down in the presence of mo- 
bile charges. 
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Figure 1: (a-c) Disorder-induced real-space magnetization from a single ran- 
dom weak impurity configuration with B = 0, V"" 1 ' = 3t, T = 0.025?, 
n imp = x = 7.5% and U = 2.8f (a), U = 3.2t (b), and U = 3.6t (c). (d) 
Schematic illustration of how a single Neel phase gets stabilized by disorder. 



The present mean field treatment of the Hubbard interac- 
tion does not capture the band narrowing due to correlations 
near half-filling leading to the Mott transition. However, it is 
expected that underdoped systems are characterized by larger 
effective interactions. Within the present scenario for dirty 
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cuprates, the x dependence of the spin order is therefore de- 
scribed qualitatively by the sequenceQJc) — >|TJb) — >QI a X until 
it disappears completely at effective U's below U c \ near op- 
timal doping. Increasing x should also be accompanied by a 
weakening of the dopant disorder potential due to enhanced 
screening. Within our model, increasing U or V" np leads to 
qualitatively similar results, and we cannot determine from this 
approach which effect is dominant in real systems. Note from 
Fig. [T]that AF droplets are induced for U > 2 At, a substantially 
reduced critical value compared to the 1 -impurity case where 
U c \ = 3.25f. This is because the Hubbard correlations induce 
charge redistributions which alter the effective local chemical 
potential, such that the criterion for magnetization of each im- 
purity site depends on its local disorder environment. Some re- 
gions containing impurities have charge densities closer to the 
phase boundary for AF order, thus enhancing local moment for- 
mation relative to the single impurity case. In the limit of large 
U, the magnetic order becomes qualitatively similar to that aris- 
ing in a stripe state with quenched disorder[40, 54, 55, 56]. 

In the presence of a magnetic field, vortices may lead to an 
enhancement of the SDW order. This is illustrated in Fig. [2] 
which shows the magnetization for a random disorder poten- 
tial with 5 = Eh)] and B + Etc)]. In Fig. Etc) two flux 
quanta penetrate the system corresponding to a large magnetic 
field of B — AIT . Focussing on the T-dependence, Fig. |3a) 
shows the structure factor \M(q)\ 2 integrated near q = (n,n) as 
a function of -T|Q]]. The main result of Fig. [2f a) is seen by the 
dashed and dot-dashed lines displaying \M(q)\ 2 including disor- 
der with and without an applied magnetic field, respectively. In 
both cases the magnetic order sets in at T c below which the T- 
dependence is in remarkable agreement with neutron diffraction 
data on underdoped LSCO|Ql. The origin of the qualitatively 
different curvature of the T-dependence of \M(q)\ 2 in the pres- 
ence of a magnetic field can be traced to a larger concentration 
of anti-phase domain walls when B = 0|uJ]. 

We now turn to a discussion of the transverse dynamical sus- 
ceptibility x^(fi, ?), oS) = -i^dte^ ([cr c (?),cr c (0)]) which 
can be expressed in terms of the BdG eigenvalues E„ and eigen- 
vectors u„, v„ 



X x x (?i,?j,aj) 
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where f(u, v) = u* mJ v* nJ [u m jV„j - u„ 7 jV mj j, and g(u, v) = 

v m iu„ i (u* ; v* . - u* v* ). Including the electronic interactions 
within RPA we find for the full susceptibility 
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Fourier transforming with respect to the relative coordinate 
f = ?i — fj defines the spatially resolved momentum-dependent 

susceptibility x(q, R, co) - e ' c ' r x(R> w )- Averaging over the 
center of mass coordinate 3 — (r, ■ + rj)/2, this expression gives 
the susceptibility x(q, co) relevant for comparison with neutron 
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Figure 2: (a) -"-dependence of the peak intensity of the magnetic structure fac- 
tor integrated around (n, n) in zero field and finite field with the density of im- 



purities rij, 



10% (dash-dotted and dashed curves, respectively). For 



the data with d> = 2<_>o and n„ 



the zero-field data were subtracted. For 



comparison, the structure factor in a clean system is also included for the same 
magnetic-field strength (solid curve). \M(q)\ 2 (integrated) translates directly to 
the ordered spin moment squared in units of fig per Cu 2+ ion. (b,c) Real-space 
magnetization at T = 0.025f for a system with 10% weak impurities modeling a 
dopant potential (b) and with the addition of two vortices (c). In (c) the vortices 
are located at sites (7, 6) and (15, 17). For all results in this figure U = 2.9t. 



scattering measurements. We have checked that at half-filling 
in the clean AF state without pairing, xiq, co) is characterized by 
spin-wave excitations on a cone centered at (n, n) as expected. 

Although we are restricted to a mean-field/RPA approach, 
a strength of the present calculation is that we can include 
both spin and charge degrees of freedom within an unrestricted 
method which is capable of describing realistic inhomogeneous 
situations, and thus naturally includes both the response from 
regions dominated by the dSC condensate or local SDW order. 
For a homogeneous dSC our model reduces to a system of Bo- 
goliubov quasiparticles whose magnetic response displayed in 
Fig. [3ja) depends crucially on the presence of a dSC gap. In the 
remainder, we investigate the role of spatially inhomogeneous 
local moments induced in this quasiparticle system. Figure[3jb) 
shows how the magnetic droplet phase freezes the low-energy 
spin fluctuations by removing the spin-gap of the hour glass 
spectrum[40]. This is also seen from the ^-summed local sus- 
ceptibility xiw) - IZqXiQ'^) plotted in Fig. [3jc) for different 
concentrations of the nonmagnetic disorder. In Fig. [3jb) and 
[3]c) we have not included the weaker disorder potential aris- 
ing from the dopants which will lead to additional freezing, but 
only included disorder from a few percent strong scatterers rel- 
evant for the modeling of e.g. Zn substitutional disorder. In 
Figs. Oa,b,c) the spectrum is C\ symmetric both in the static 
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and dynamical response. In Ref. H40fl is was shown how the 
stripe phase breaks this C\ symmetry in the inelastic response 
even in the presence of substantial disorder concentrations. 




co/t wit 



Figure 3: (a,b) Spin susceptibility x(l, t, <u) for a clean (dirty) dSC phase with- 
out (with) disorder using U = 2. (it, B = 0, T = 0.04/, and in (b) n imp = 3% 
strong scatterers with V""'' = lOOr. (c) Local susceptibility x^) f° r the same 
parameters as in (b) except for different disorder concentrations of strong scat- 
terers. In (b,c) we have averaged over 5 different impurity configurations, (d) 
comparison of xt 10 ) in the clean dSC phase with (without) an applied magnetic 
field shown by dashed (solid) line. 

We end by discussing the dynamical spin susceptibility when 
B + 0. Since the vortices can also nucleate local SDW order, 
we expect a similar slowing down of the fluctuations as for the 
point-like disorder. That this is the case can be seen from Fig. 
[3d) which compares ^(w) for the clean dSC with and without 
an applied magnetic field[57]. Clearly, the same filling-in of the 
low-energy weight takes place because of the vortices. 

4. Conclusions 

We have reviewed a theoretical scenario for the understand- 
ing of disorder- and magnetic field-induced SDW phases in un- 
derdoped cuprate superconductors. In the weak coupling limit, 
nonmagnetic impurities nucleate static SDW order. The li-wave 
nature of the pairing state is crucial for generation of local mag- 
netism and works as a catalyst for the SDW phase. In the strong 
coupling limit disorder mainly acts to scramble correlation- 
generated stripe order. We demonstrated explicitly how both 
disorder and magnetic fields slow down the spin excitations as 
observed in experiments. 
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